Scaling approach to order-parameter fluctuations in disordered frustrated systems 
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We present a constructive approach to obtain information about tlie compactness and shape of 
large-scale lowest excitations in disordered systems by studying order-parameter fluctuations (OPF) 
at low temperatures. We show that the parameter G which measures OPF is 1/3 at T = provided 
the ground state is unique and the probability distribution for the lowest excitations is gapless and 
with finite weight at zero-excitation energy. We then apply zero-temperature scaling to describe the 
energy and volume spectra of the lowest large-scale excitations which scale with the system size and 
have a weight at zero energy Pv(0) ~ with v = l'^. A low-temperature expansion reveals that, 
OPF vanish like , ii 9 > and remain finite for space filling lowest excitations with ^ = 0. The 
method can be extended to extract information about the shape and fractal surface of the large-scale 
lowest excitations. 



Introduction. There is a long standing controversy on 
the nature of low-temperature excitations in spin glasses. 
This discussion has seen a strong revival in the last years 
due to the large increase of computational capabilities 
together with the simultaneous developnjiei|it of more re- 
fined and powerful numerical techniquesun. These have 
largely improved the ability to find exact ground states 
or to equilibrate small systems at very low temperatures. 
The aim of this letter is to present some analytical results 
on this subject by investigating order parameter fluctua- 
tions (hereafter denoted by OPF) and showing how they 
can give us information about the topology of the exci- 
tations in spin glasses. 

The study of parameters measuring OPF has been 
proven to be very powerful to locate fiflte temperature 
phase transitions in disordered systems Recent stud- 
ies □ show that these parameters can bring also informa- 
tion about the frozen phase and thus shed some light on 
the nature of the spin-glass phase. The introduction of 
G was originally motivated by the existence of the well- 
known Guerra sum rules for spin-glassea3, and in partic- 
ular the following one: (q^)^ = ^(9^) + |('?^) where 
q ^ y Si^i '^i'^i is the usual overlap between different 
replicas {c, r}, (...) and (...) stand for thermal and dis- 
order averages respectively. This relation implies that in 
systems with non vanishing OPF in the frozen phase, the 
parameters G and A defined as 



(1) 



have the following behavior in the infinite- volume limit: 
G = l9{Tc - T) and A = A{T) is a non-vanishing func- 
tion of T (this result for A being an evidence of replica- 
symmetry breaking (RSB)). ^ 

Very recently it was conjectured in Q (henceforth re- 
ferred to as RS) that for systems of Ising variables G 
takes the universal value 1/3 at T = for finite volume 
systems under two mild assumptions: the uniqueness of 
the ground state and the absence of zero gap in the local 



field distribution. In RS the main scenarios for OPF were 
discussed and many examples were given to check the va- 
lidity of such a conjecture. Moreover, the conjecture was 
proven for finite volume systems provided that the limit 
T — > is taken before the limit F ^ 00 by showing that 
one-spin excitations with finite probability at zero local 
field yield at low temperature the dominant (linear in T) 
contribution to OPF. 

In this letter we prove that these results generally hold 
whatever the size of the lowest excitations is, provided 
these are gapless and have finite weight at zero excita- 
tion energy. Moreover we will obtain valuable informa- 
tion about the size of the excitations by considering the 
low-temperature behavior of the OPF parameters. In 
particular, we will see how we can establish a scaling 
theory for the lowest large-scale excitations which relates 
OPF to the T — Q thermal scaling exponent 0' and the 
exponent which characterizes the distribution of large- 
volume excitations 7. We exemplify these findings in the 
spin-glass chain. In this paper we analyze the excitations 
through the overlap parameter q and we do not discuss 
about the shape of the excitations, an issue which be ad- 
dressed by studying the energy overlap or link overlap 
order parameters in a similar way. 

The model. Consider a system of V Ising spins with a 
general Hamiltonian 



n = - ^Jij(Ji(Tj 



(2) 



where the J^j are the quenched variables. In general 
these can be either (depending on the dimensionality 
and interaction range of the system) or distributed ac- 
cording to some continuous P{J). For a continuous dis- 
tribution of couplings there is a unique ground state (up 
to a global flip of all the spins), whose configuration we 
denote by {cr*}. This Hamiltonian can be rewritten in 
terms of local fields, following the formulation given in 
RS, Ti. — ~ CTj hi . In the ground state configuration 
each spin ai is aligned with its local field hi, so that the 
ground-state energy is £:gs = ~ Ei = - Ei 1^*1 = 
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" Si j '^ij'^i'^j- Suppose that for a given sample we know 
the ground state configuration {a*}. In this state, we 
will take the convention that all the spins have the same 
color, let say red. Consider now a cluster excitation of 
volume V. In this new state v spins will change their ori- 
entation with respect to the ground state configuration 
so that they have a different color, let say blue. Thus, 
in the excited state the spins will fall into two groups: 
the red one, which we denote by V, and the blue one, 
which we denote by V. This minimum cost excitation is 
always a cluster excitation in the sense that all excited 
spins are connected to each other by bonds. An excita- 
tion consisting in two or more disconnected clusters will 
necessarily cost more energy than reverting one single 
cluster. Since we are considering a cluster excitation, we 
can define among these two regions the surface of the ex- 
citation or droplet: 5V. SV will be constituted by those 
pairs of spins of different color which are connected to 
each other by at least one bond. Thus the surface con- 
tains red ((5P fl V) and blue spins {SV D V). 

The energy cost of such an excitation will depend ex- 
clusively on the bonds which connect spins of different 
colors, defining the surface of the droplet. Therefore, the 
energy cost reads, 

A£; = 2 J^^k<J:al , (3) 

i,k diff. color 

{i,k)eSV 

where each coupling is counted only once. It follows that 
the lowest excitation with energy {AE*) will correspond 
to breaking a subset of bonds among the {Ji,k} such that 
the value in (|^) is minimum. 

Analysis of OFF fluctuations. If we are at a very low 
temperature the main contribution to the partition func- 
tion arises from the lowest excitation. In this limit the 
partition function reads 

The computation of the terms involved in G is done 
following the same procedure used in RS. Different quan- 
tities which enter in the definition of G can be expressed 
in terms of the two and four point correlation func- 
tions by defining the following objects = (aiOj)^ and 

Fijkl = {(^if^jCTkOl)'^ ■ 

i^^j (ij,k,l) 

To compute the two-point correlation functions {(JiOj) 
we must take into account two different contributions 
which we denote as Uij and Vij: 

• Pairs of spins of the same color, either blue tTj^) = 
— cr*, or red (Ji^j) = cr* , both yield the same contri- 
bution {(JiOj) — cr*a* = Uij. 



• Pairs of spins having different colors yield [aiOj) = 
a,ra;tanh(^) ^1^,, . 

In the same way the four-point correlation functions 
{(JiajakCTi) can be computed by adding two types of con- 
tributions: Uijki and Vijki- 

• When either, all four spins have the same color 
or there are two spin pairs with different colors, 

• When three spins have the same color and the 
other spin has a different color {aiCFjakcri) = 

<^*3^k^i tanh = V^ju . 



The next step in order to compute {q^), (<z^)^ and (q"*) 
is to take into account the degeneration dy of each kind of 
contribution. This number will depend on the size of the 
excitation (w). Consider, for instance, the term X^i/j -^u - 
An excitation with v blue spins and {V — v) red ones con- 
tributes to Tij with the following degenerations: 

• The number of pairs of spins which in the excited 
state are both red {d^ — {V — v){V — v — 1)) or 
both blue [d^ = v{v - I)) yield Uf^ = 1. 

• The number of pairs of spins with different colors 
(d„ = 2{V ~ v)v) yield V^]=R = t&nh^{^^). 

Adding all contributions we finally obtain: X^i^i^j -^u ~ 
{V - v){V -v-l) + v{v-l) + 2{V - v)vR . Up to now 
we have contented ourselves with computing the corre- 
lation functions but we have not made any hypothesis 
on the probability distribution of having an excitation of 
certain volume (u) and energy cost (AE'*). When per- 
forming the average over the disorder, we have to bear 
in mind that each sample may have a lowest excitation 
of size V and energy AE* , yielding a distribution of the 
lowest excitation energies for every size v. Since there 
may exist correlations among fields which enhance clus- 
ters of certain sizes we make the assumption that ex- 
citations of volume v {< V) can take place with a fi- 
nite probability 9v — 1) and the energy cost of 
these excitations at fixed volume v is distributed accord- 
ing to some normalized continuous Py{AE*), which is 
only defined for positive AE* . Through these considera- 
tions we get the following general expression X^i^^j — 

Y.v=i 9v [V{V - 1) + 2v{v -V) + 2{V - v)v R^] where 

i?^=/dAS*n(A£;*) (tanh{^^))"' . 

In a similar way we can evaluate all terms entering in 
. After some algebra we obtain for the numerator (N) 
and denominator (D) of G the following expressions. 
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V 

V 

V 

-{Y,g,(Rl^l)v{V-v)Y] . (7) 

V 

Note that these expressions are invariant under the 
change v ^ V — v as it should be (for every excitation 
there is the corresponding one related by time-reversal 
symmetry). Now we compute the low-temperature limit 
of the quantities i?J with k — 1,2. A simple calcula- 
tion yields, in the limit T — > keeping terms up to 
first order in T, ^ = 1 - 2Tpy{0) + 0{T^) = 
1 - |rP„(0) + 0{T^) . Substituting them into expres- 
sions (^, (^, we obtain that both numerator and de- 
nominator vanish linearly with T but its ratio remains 
finite yielding the expected result, 

G=\ + 0{T) , (8) 
A=^,Y.9M0)^'{V-vf + O{T^) . (9) 

v — l 

The only necessary condition to obtain these last re- 
sults is that Py (0) has a finite weight at zero energy cost 
for one given volume (u). Otherwise, the expression for 
would have corrections of order 0{T^) and the uni- 
versal value 1/3 for G and the lineai^in T dependence of 
A would not be recovered anymordZl. Let us note that 
in RS we considered one-spin excitations so we assumed 
that, P^(0) ^ for v = 1. In that work we stressed that 
in order to have corrections 0(T) when considering the 
inversion of two or more spins there had to be singular 
correlations between local fields at different sites. This is 
the situation we expect to meet in generic spin-glass sys- 
tems where large cluster excitations are responsible for 
the universal low-temperature value of G. 

Scaling theory for OFF fluctuations. Up to now we 
only considered the linear terms in a low T expansion. 
But this first order term suffices to establish the scal- 
ing behavior of A. In this approach we will focus on 
large-scale excitations which scale with the volume of the 
system, assuming that they bring the relevant contribu- 
tion to thermodynamic quantities. Let's suppose that 
the weight at zero energy of lowest excitations scales 
like Py^id{0) ~ , where d is the dimensionality of 
the system and 9' is the thermal exponent for the low- 
est excitation^!. Then we suppose the following scaling 
ansatz for lowest-lying large-scale excitations: 

g,^L-^g{^) ,P„^;.(0)^r«' ,-^OC0(l) . (10) 

where L is the lattice size and c? < 7 (since g^ is posi- 
tive defined and normalizable) . This functional form im- 



plies that in order to assure the normalization condition 
J^vdi'") ~ weight of finite volume excitations is 

1 — L~'^. Note that, in principle, since excitations of any 
size are possible, the average volume of the excitations 
scales like v '--^ L'^'^~'' 113. Substituting these relations in 
A in (H) and taking the large L limit we get, 

A = ^TL-^+'^-«' dxg{x)x^-'^ (l + . (11) 

Assuming that g{x) decays fast enough for a; ;2> 1 (e.g. 
exponentially), the general low T scaling A ^ TL~^, 
6 = 9' — d + ^ together with, relations ( p^ ) provide a way 
to estimate both 9' and tIIII. In this scenario, since A is 
bounded A < 1, the crossover behavior of A is then de- 
termined by the term TL~^ which gives a typical length, 
scale I* ~ T^^^. In usual domain wall scaling theoryS 
we expect that, for 9 < 0, there is a relation between 
the thermal exponent 9 which accounts for the scaling of 
the free-energy of the largest excitationaa, and the cor- 
relation length exponent v = —1/9. Below the lower 
critical dimension di we have 6* < 0. Above di, 9 > and 
the low-temperature behavior of A may decide between 
the two most rcontroversial spin-glass scenarios. In the 
droplet modeB 6* > and A - T/L^ while in the RSB 
scenario 9 = yielding a finite A at low T in the L ^ 00 
limit. This result corresponds to the case where there is a 
finite probability density of samples with large scale exci- 
tations of finite energy implying that Py{AE*) for v ^ V 
is scale independent. The scaling relation A ~ T/L^ can 
be checked by numerical simulations through measure- 
ments of A (preliminary pjsuJts in 2 and 3 dimensions 
agree with that expressior£3ll3) . We remind that it holds 
only for T/L^ ^ 1 (where the approximation (^ is valid) 
and does not invalidate the present analysis for 9 < ii 
T is low enough. 

The Ising spin-glass chain. An illustratiszje example of 
these results is the Ising spin-glass chairo which is de- 
scribed by the following Hamiltonian: 

L 

H = —'^^JiCJiUi+i , (12) 

where the quenched disorder is distributed according to 
some continuous distribution P{ J) with finite weight at 
zero coupling. The ground state is unique (up to a global 
spin inversion) and has an energy Eqs = — J2i Ji<^i(^i+i- 
In the study of the low-temperature excitations we have 
to distinguish between the cases with free and periodic 
boundary conditions. In the free boundary case the mini- 
mum cost excitations correspond to breaking the weakest 
bond J* and reversing all the spins on the right or left of 
that bond. The energy cost of such an excitation reads: 
IS.E* = 2\J*\. Note that depending on the sample, the 
length of the droplet (w) can range between 1 an L — 1. 
Clearly these excitations have a finite probability of hav- 
ing zero cost as P{l\E* = 0) - L p\j* = 0) which is 
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finite by hypothesis. Hence from our previous general 
arguments it follows that at T = we obtain the univer- 
sal value G = 1/3 in agreement with the analytical result 
derived in RS. It is important to realize that one-spin ex- 
citations cannot correspond to minimum cost excitations 
when considering spins which are not at the boundaries. 
This is due to the fact that, in this particular system, 
the ground state is not frustrated. Exciting one interior 
spin would necessarily frustrate two bonds so that the 
energy cost would be: Ai?* = 2(|Ji| -I- |Jfc|) , which is 
obviously less favorable than breaking a single bond. It 
is easy to check that Sv = substituted in (^) yields 

A = in full agreement with transfer matrix 

calculations shown in RS. Eqs. ( [lO| , [ll| ) with 7 = 0?=! 
yield 9' = -1. 

In the periodic boundary conditions case, the situation 
is different. Here we have L spins and L couplings, 
with the condition ctl+i = ci- Note that the ground 
state can have frustrated bonds, so that we can dis- 
tinguish between two different kinds of samples: non- 
frustrated (NF) having Yii Ji > and frustrated (F) 
having Yii Ji < 0. In the first case the ground-state 
energy reads: E^p = \Ji\ while in the second we 

have that the weakest bond, let's say Jk, is unsatisfied so 
that the ground-state energy reads: Ep = E^p + 2| J^j. 
In both cases, the minimum energy excitation will corre- 
spond to breaking the two weakest bonds, let's say (z, fc), 
but the energy cost will read differently in each case: 
A^i^F = 2(1 J,| + iJfcl), A£;^ = 2(1 J,| - iJfcl). From here, 
it is evident that the energy gap in the F set will be 
much smaller than in the NF set. Thus, only F_aam- 
ples yield contributions 0{T) at low temperatures Il3. In 
other words, P{AE§. =0)^0 while P{AE;^p = 0) = 
ahhough, in average, still G(T = 0) = 1/3 and A TL. 
We conclude summarizing our results. A low T ex- 
pansion for OFF reveals that the only condition to get 
G = 1/3 + 0{T) and A - T is that the ground state 
is unique and that the disorder average probability dis- 
tribution for the lowest excitations is gapless and with 
finite weight at zero-excitation energy. This illustrates 
the importance of rare samples in determining the low- 
temperature properties of spin glasses. Assuming that 
the statistics of the large-scale lowest excitations deter- 
mines the low T thermodynamic behavior usual scal- 
ing arguments yield -?V=L''(0) ~ and A ^ TL^^ 
providing a way to estimate 9' and 9. In this direc- 
tion, numerical investigation of the . firat lowest excita- 
tions in two-dimensional spin glassest§ll3 would be wel- 
come to confirm the validity of the scenario we have put 
forward through the scaling relations (|o|). Let us also 
mention that most of the predictions presented here (in 
particular relations (^,^,|ll])) can be checked with nu- 
merical simuiations for small sizes thermalized at low- 
temperatureaa. The calculations shown here open other 
possibilities. Here we concentrated our attention in the 
study of the parameter q which gives information about 
the size of the excitations. To investigate other topolog- 
ical aspects of the excitations (for instance, their fractal 



surface exponent ds) a similar study of OFF for energy or 
link overlaps along the present lines should be necessary. 
Further information about the topology of the excitations 
can be obtained by imposing G = 1/3 to all orders in T, a 
result which apparently holds in spin glasses in the large 
volume limit. 
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